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ABSTRACT 

Recent cosmological observations indicate that the reionized universe may have started at around 
z = 6, where a significant suppression around Lya has been observed from the neutral intergalactic 
medium. The associated neutral hydrogen column density is expected to exceed 10^^ cm^^, where it is 
very important to use the accurate scattering cross section known as the Kramers-Heisenberg formula 
that is obtained from the fully quantum mechanical time-dependent second order perturbation theory. 
We present the Kramers-Heisenberg formula and compare it with the formula introduced in a heuristic 
way by Peebles (1993) treating the hydrogen atom as a two-level atom, from which we find a devitaion by 
a factor of two in the red wing region far from the line center. Adopting simple cosmological models, we 
compute the Gunn-Peterson optical depths and the trough profiles. Our results are compared with the 
works performed by Madau & Rees (2000), who adopted the cross section introduced by Peebles (1993). 
We find deviations up to 5 per cent in the Gunn-Peterson transmission coefficient for an accelerated 
expanding universe in the red off-resonance wing part with the rest wavelength AA ~ 10 A. 

Subject headings: cosmology: theory — intergalactic medium - quasars: absorption lines — radiative 
transfer — galaxies: high-redshift 



1. INTRODUCTION 

The quasar absorption systems have been excellent tools 
to investigate the intergalactic medium (IGM), from which 
it has been well known that the IGM of the nearby uni- 
verse is highly ionized (e.g. Peebles 1993). Since the 
universe after the recombination era z ~ 1100 should be 
dominantly neutral, there must be some epoch when the 
universe began to be re-ionized. Intensive studies have 
been performed on the emergence of the first objects that 
ended the dark age of the universe. Numerical calcula- 
tions adopting the cold dark matter models predicted the 
reionization epoch in z ~ 6 — 12 (e.g. Gnedin & Ostriker 
1997). 

Around this epoch a broad absorption trough in the 
blue part of Lya is expected and regarded as a strong 
indicator of the reionization of the universe, which was 
predicted by Gunn & Peterson (1965) and independently 
also by Scheuer (1965). With the advent of the Hubble 
Space Telescope and 8 meter class telescopes there have 
been extensive searches for the Gunn-Peterson trough in 
the spectra of high red shift objects. Remarkable contri- 
butions are made by the Sloan Digital Sky Survey, from 
which a number of high red shift quasars with z ranging 
from 4 to 6 have been found. According to the recent re- 
port from the Keck spectroscopy of these high red shift 
quasars (Becker et al. 2001), the flux level drop around 
Lya is much higher for the quasar with z — 6.28 than 
those for other quasars with z < 6, which indicates that 
the reionization epoch may be found at around z ^ 6. 

The exact computation of the flux drop around Lya re- 
quires an accurate atomic physical estimation of the scat- 



tering cross section. Recent theoretical works on the cal- 
culation of the Gunn-Peterson trough were provided by 
Miralda-Escude (1998) and Madau & Rees (2000), who 
adopted the formula that was introduced in a heuristic 
way using the second order time-dependent perturbation 
theory by Peebles (1993). The formula is derived based 
on the assumption that the hydrogen atom is a two-level 
atom, in order to show the behavior of the scattering cross 
section that is approximated by the Lorentzian near res- 
onance and yields dependence in the low energy limit. 
As Peebles noted clearly in the text, due to the two- level 
assumption, the formula provides an inaccurate propor- 
tionality constant in the low energy limit, even though it 
correctly gives the dependence. 

The accurate cross section should be obtained from the 
second order time-dependent perturbation theory treat- 
ing the hydrogen atom as an infinitely many-level atom 
including the continuum free states, which is known as 
the Kramers-Heisenberg formula. The discrepancy be- 
tween the two formulae will be significant in the far off- 
resonance regions where the contribution from the np{n > 
2) states including the continuum states becomes consid- 
erable. Therefore, in order to obtain an accurate Gunn- 
Peterson profile it is essential to investigate the exact scat- 
tering optical depth of a medium with a high neutral hy- 
drogen column density. 

In this Letter, we present a faithful atomic physics that 
governs the scattering around Lya by introducing the 
Kramers-Heisenberg formula and a simple fitting formula 
around Lya. We compute the Gunn-Peterson trough pro- 
files adopting a representative set of cosmological param- 
eters with our choice of the reionization epoch and make 
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quantitative comparisons with previous works. 

2. THE KRAMERS-HEISENBERG FORMULA 

The interaction of photons and electrons is described by 
the second order time-dependent perturbation theory, of 
which the result is summarized as the famous Kramcrs- 
Heisenberg formula. As is well illustrated in a typical 
quantum mechanics text, it may be written as 
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where e" , e" arc the polarization vectors associated with 
the incident photon and outgoing photon, respectively, F/ 
is the radiation damping term associated with the inter- 
mediate state 7, w is the incident angular frequency, is 
the angular frequency of the transition between I and the 
ground state A, and Tq = e^/{meC^) = 2.82 x 10^^'^ cm is 
the classical electron radius (e.g. Sakurai 1967). 

Here, the electron is in the ground state A before scat- 
tering and de-excites to the same state. The summation 
(and integration) should be carried over all the intermedi- 
ate states / including the infinite number of bound states 
and the free or continuum states. For hydrogen, the dipole 
moment matrix elements have been explicitly given using 
the reciirrence relations of the hypergeometric func;tion in 
many texts and in the literature (e.g. Berestetski, Lifshitz 
& Pitaevskii 1971, Bethe & Salpeter 1957, Karzas & Lat- 
ter 1961). 

In the blue part of Lya, it is possible that the scattering 
atom may de-excite to the excited 2s state by re-emitting a 
photon with nnich lower frequency than the incident pho- 
ton. This inelastic scattering or the Raman scattering is 
negligible near Lya due to small phase space available for 
an outgoing photon. However, this process becomes im- 
portant as the incident photon energy increases. In the 
range where the current work is concerned, the Raman 
scattering process is safely neglected. 

Despite the existence of the explicit analytic expres- 
sions of each matrix element that constitutes the Kramers- 
Heisenberg formula for hydrogen, it is still cumbersome to 
use the formula as it is. Therefore, a simple fitting for- 
mula around Lyo will be useful for practical applications. 
Near resonance (1170 A < A < 1410 A), the Lorentzian 
function gives quite a good approximation 
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where F2p = 6.25 x 10^ s~^ is the radiation damping con- 
stant associated with the Lya transition. In the long wave- 
length region (A > 1410 A), Gavrila (1967) provided the 
fitting polynomial for the Rayleigh scattering cross section, 
which is 

a{uj)/aT = QAm{Lo/uJaf + 0.900(a;/a;a)^ 

+ \2.&{uj/uJo.f\ (3) 

In particular, Ferland (2001) applied Gavrila's fit to his 
photoionization code 'Cloudy'. 



In the case of the short wavelength region (1070 A < 
A < 1170 A), we provide a similar polynomial fit to the 

Kramers-Heisenberg formula 

cr(w)/c7T = 1.62 X l(f{uja/iof + 5.88 x l{f{uja/ujf 
+ 7.99 X 10^(w«/w)^ - 4.83 x lQ^{uj„/w) 
+ 1.09 X 10'^. (4) 

The deviation of the fit is within 5 per cent from the true 

Kramers-Heisenberg formula. 

In Fig. 1 we show the scattering cross section from the 
Kramers-Heisenberg formula by the solid line and by the 
dotted line we represent the fit. The behavior near reso- 
nance is depicted in the bottom panel, because the cross 
section changes very steeply. It is apparent that the scat- 
tering cross section is excellently approximated by the 
Lorentzian. However, in the wavelength range considered 
in this Letter, the radiation damping is completely neg- 
ligible, and the curve shown in the figure is simply pro- 
portional to Ao)^^ = (cl) — (jJLycx)^"^- The deviation is 
slightly anti-symmetric with respect to the line center in 
the sense that the cross section in the blue part is smaller 
than the Lorentzian and in the red part it is larger than 
the Lorentzian. 

Therefore, 1 per cent of deviation of the Lorentzian from 
the Kramers-Heisenberg formula is seen at a wavelength 
shift of AA = A — Aq, = ±3.3 A, for which the correspond- 
ing cross section a = 3.8 x 10^^^ cm^. This indicates the 
accuracy of the Voigt profile fitting applied to quasar ab- 
sorption systems, where the accuracy is more than 99 per 
cent when the absorbing medium is characterized with the 
H I column density smaller than 3 x 10^*^ cm~^. 

Further away from the line center, the cross section in 
the blue part decreases very steeply till A ~ 1100 A, but 
the decrease of the cross section in the red part is rather 
gradual and eventually becomes proportional to w^, which 
corresponds to the classical result. 

Fig. 1 also shows the comparison of the Kramsers- 
Heisenberg formula and the heuristic formula 
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introduced by Peebles (1993). 

The dependence in the limit a; <C is obtained as a 
result of the closure relation, which is apparent in the both 
formulae. However, the Kramers-Heisenberg formula gives 
about twice larger scattering cross section than ap does. 
This deviation is easily noted when the oscillator strength 
of the Lya transition is /lj/q = 0.42. The scattering cross 
section in the far red region is contributed from all the p 
states and the oscillator strength is a good measure of the 
contributions of each individual excited state. This im- 
plies that the 2p state contribution is comparable to the 
total contributions from the remaining states in the low 
energy limit. 

3. GUNN-PETERSON TROUGH PROFILES WITH THE 
KRAMERS-HEISENBERG FORMULA 



We compute the Gunn-Peterson optical depth defined 



by 



TGP = I dz ^^[v = c(l + z)/\obs] n{z), (6) 
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with the Kramcrs-Hciscnberg formula and make compar- 
isons with previous works performed by Madau & Rees 
(2000) (see also Miralda-Escude 1998). Here, Xobs is the 
observed wavelength, Zrd, Zg are the redshifts of the com- 
plete reionization of the universe and the rcionizing source, 
and n{z) = no{l + z)^ is the homogeneous neutral hydro- 
gen density at redshift z. We choose Zrei = 6, = 7 as in 
Madau & Rees (2000), but do not consider the proximity 
effect of the ionizing source. The Gunn-Peterson optical 
depth can be written as 

Tap = NoHi dz a[v = c(l + z)/\obs] 

x{i + zf/[nM{i + zf + niyY'\ (7) 

where Q,m,^k are the density parameters due to matter 
and the cosmological constant and the characteristic hy- 
drogen column density 

NoHi=nocHoK (8) 

Wc choose the present Hubble constant and the hydro- 
gen number density Hq = 50 km s^^ Mpc~^, no = 2.4 x 
lO"'^ cm-^ so that Nqhi = 4.3 x 10^^ cm'^. 

In Fig. 2, we show the Gunn-Peterson transmission 
coefficient Tqp = 6""^°^. We plot Tap for the case 
Hm = 1, = in the top panel and the same quantity for 
an accelerated expanding universe ^Im = 0.35, flA = 0.65 
in the bottom panel. In terms of the characteristic Gunn- 
Peterson optical depth tqgp defined as 

3Xlj:2pn{zs) 
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our choice of parameters in the case of the top panel yields 
the value of tqgp = 3 x 10^ at Zg = 7 as was adopted in 
the work of Madau & Rees (2000). 

The horizontal axis represents the logarithm of the nor- 
malized wavelength ratio S defined as 

A,,,,, 



from which 6 = corresponds to the resonance wavelength 
of Lya. By the dotted lines we present the Gunn-Peterson 
transmission coefficient obtained using the scattering cross 
section given by Eq. (5). The deviation between the two 
formulae is notable around S = 10~^, where the deviation 
is about 3 per cent in the top panel and 5 per cent in the 
bottom panel. The Peebles approximation turns out to be 
pretty good for contemporary application. 



The deviation between the two formulae will increase 
as no or Nqhi increases, because the discrepancy of the 
Kramers-Heisenberg formula and Eq. (5) becomes larger 
as the frequency is further away from the line center. Near 
resonance, both formulae are excellently approximated by 
the same Lorentzian. Therefore, no significant deviation 
is expected when the neutral medium is of low column 
density < 10^^ cm~^. It is notable that an accurate treat- 
ment of atomic physics is more important in an accelerated 
expanding universe where the univere was more compact 
than the universe without the cosmological constant. 

4. SUMMARY AND DISCUSSION 

In this Letter, we have investigated the behavior of the 
scattering cross section around Lya in a quantitative way, 
where the deviation from the Lorentzian becomes signifi- 
cant as the incident frequency gets further away from the 
line center. Therefore, in an analysis of the Gunn-Peterson 
trough profile, which is associated with a neutral medium 
with a high H I column density, an inaccurate treatment of 
the atomic physics of hydrogen may introduce significant 
errors in estimating important cosmological paramet(>rs in- 
cluding the epochs of the emergence of the first objects and 
the completion of the reionization of the universe. 

Voigt profile fitting has been very successfully applied 
to quasar absorption systems with a broad range of H I 
column densities. However, the deviation of the true scat- 
tering cross section from the Lorentzian exceeds 1 per 
cent when the relevant column density becomes Nhi > 
3 X 10^" cm~^ that is the typical column density of a 
damped Lya absorber. This is especially important in 
some damped Lya systems that may possess Njji > 
10^1 cm-2 (e.g. Turnshek & Rao 1998). However, it 
should be noted that the damping constant T2p is so small 
compared with the scale relevant in this work, the cross 
section is effectively of the form oc Aw"^ = (w — Wa)"^. 
Therefore, the absorption profile is irrelevant to the exact 
value of the radiation damping term, which means that 
the term 'damped Lya absorption' is a misnomer. 



KHB thanks the support from the BK21 project initi- 
ated by the Ministry of Education. HWL thanks Roger 
Blandford and Sang-Hyeon Ahn for their helpful discus- 
sions. 



REFERENCES 



Becker, R. H., Fan, X., White, R. L., Strauss, M. A., Narayanan, V. 

K., 2001, AJ submitted (astro-pli/0108097) 
Berestetskii, V. B., Lifshitz, E. M., & Pitaevskii, L. P., 1971, 

Relativistic Quantum Theory, Pergamon Press, Oxford 
Bethe, H. A., Salpeter, E. E., 1957, Quantum Mechanics of One and 

Two Electron Atoms, Academic Press Inc. New York 
Fcrland, G., Hazy, a brief introduction to Cloudy 94-00 
Gavrila, M., 1967, Phys. Rev., 163, 147 
Gncdin, N. Y., & Ostrikcr, J. P., 1997, ApJ, 486, 581 
Gunn, J. E., Peterson, B. A., 1965, ApJ, 142, 1633 



Karzas, W. J., & Latter, R., 1961, ApJS, 6, 167 
Madau, P., Rees, M. J., 2000, ApJ, 542, L69 
Miralda-Escude, J., 1998, ApJ, 501, 15 

Peebles, P. J. E., 1993, Principles of Physical Cosmology, Princeton 

University Press, Princeton 
Sakurai, J. J., 1967, Advanced Quantum Mechanics, Addison- Wesley 

Publishing Company, Reading, Massachusetts 
Scheuer, P. A. G. 1965, Nature, 207, 963 
Turnshek, D. A., & Rao, S., 1998, ApJ, 500, L115 



4 



Bach & Lee 




Fig. 1. — The scattering cross section around Lya. The solid Une represents the accurate, fully quantum mechanical cross section known 
as the Kramers-Hcisenbcrg formula. The dotted line represents the Lorentzian given by Eq. (2), which gives an excellent approximation 
near the line center. The long dashed line represents the fitting formula provided by Gavrila (1967), of which the approximation is valid for 
A > 1400 A. The dot-dash line represents our fit to the Kramers-Heisenberg formula in the blue part given in Eq. (4). The dot-long dash line 
represents the cross section obtained from Eq. (5), which is inaccurate by a factor of two in the far red wing region and accurate near the line 
center. In the bottom panel, we plot the cross section obtained from the Kramers-Heisenberg formula (solid line) and the Lorentzian Eq. (2). 
The cross section is asymmetric relative to the line center in the sense that the cross section in the blue part is smaller than in the red part. 
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Fig. 2. — The Gunn-Pcterson transmission coefficient Tqp = e^^f^p for r2jvf = li^A = (top panel) and for Q-m = 0.35, Qa = 0.65 
(bottom panel). The solid lines represent the values obtained using the Kramers-Heisenberg formual and the dotted lines are for the values 
from the cross section cp introduced by Peebles (1993). The present Hubble constant and the hydrogen number density are chosen to be 
ifo = 50 km s-i Mpc-i,no = 2.4 x lO"'' cm-3 so that N^hi = nocR-^ = 4.3 x lO^i cm-2. 



